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ABSTRACT 


A  comparison  is  made  between  values  of  radio  ray  bending 
which  is  calculated  in  two  ways:  one  using  a  ray-tracing  compu¬ 
tation  based  on  radiosonde  data;  and  the  other  using  the  National 
Bureau  of  Standards  exponential  reference  atmosphere  (1959) 
corrected  to  match  the  refractivity  measured  at  the  surface. 

The  refractivity  data  for  the  ray-tracing  computations  were 
obtained  from  twice-daily  radiosonde  soundings  taken  at  Albany, 
New  York,  and  Nantucket,  Massachusetts,  during  July  and 
November  1963.  Results  indicate  that  the  NBS  model  does  not 
completely  remove  height  errors,  but  it  can  be  used  effectively 
to  offset  errors  due  to  slow  seasonal  changes  in  bending  conditions. 
The  correction  technique  is  not  effective  during  periods  of  anoma¬ 
lous  propagation.  A  review  of  more  recent  prediction  techniques 
(1963)  is  presented,  and  recommendations  are  made  for  their 
field  evaluation. 


REVIEW  AND  APPROVAL 

This  technical  report  has  been  reviewed  and  is  approved. 
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SECTION  I 


INTRODUCTION 


Radio  refraction  in  the  lower  troposphere  continues  to  undergo  intensive 
investigation  due  to  the  increasing  precision  demanded  by  advanced  electro¬ 
magnetic  systems.  The  index  of  refraction  varies  significantly  in  the  vertical 
direction.  However,  Bean  and  Cahoon  have  shown  it  to  be  almost  hori¬ 
zontally  stratified  for  ray -tracing  analysis,  except  in  the  case  of  anomalous 

[2] 

conditions.  The  investigations  of  Bauer,  Mason,  and  Wilson  ‘  indicated 
that  the  vertical  variation  of  the  index  behaved  as  an  exponential  function 

[  j 

with  height  during  most  of  the  year  in  temperate  latitudes.  Bean  and  Thayer 
developed  an  exponential  reference  atmosphere  based  on  a  large  number  of 
radiosonde  soundings  from  climatically  and  geographically  diverse  U.  S. 
stations.  The  results  of  this  latter  investigation  support  the  exponential 
model  and  also  show  that  a  relationship  exists  between  the  index  at  the 
surface  and  the  average  index  gradient  measured  over  the  first  kilometer 
in  height. 


Radio  ray  bending  occurs  mostly  in  the  lower  part  of  the  troposphere 

where  the  vertical  variation  of  index  is  greatest.  Therefore,  the  fact  that 

the  initial  gradient  over  the  first  kilometer  can  be  related  to  the  readily 

obtainable  surface  value  has  great  practical  significance.  Sweezy  and 

[4] 

Bean  evaluated  the  use  of  the  surface  index,  initial  surface  gradient, 

and  first  kilometer  average  gradients  to  correct  radar  height  errors.  The 

first-order  correction  only  used  the  surface  index  N  .  The  accuracy  was 

s 

then  increased  by  using  both  Ng  and  the  initial  surface  gradient  . 
Finally,  the  analysis  included  not  only  Ng  and  Gq  ,  but  also  the  gradient 
AN/Ah  averaged  over  the  first  kilometer  in  height.  More  recently, 
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Sweezy  and  Bean  extended  their  technique  into  the  range  of  extremely 
small  initial  elevation  angles. 

Towers  and  radiosonde  soundings  would  be  required  to  determine  the 
surface  gradient  and  the  first  kilometer  average  gradient  AN/Ah. 

These  facilities  are  not  normally  or  conveniently  available  at  most  sites. 
Therefore,  this  report  investigates  the  use  of  surface  measurements  alone 
to  determine  ray-path  bending  and  height  errors  with  the  1959  NBS  model 
profile. 

Albany,  New  York,  and  Nantucket,  Massachusetts,  were  selected  for 
investigation  during  the  months  of  July  and  November  1963.  The  NBS  model 
predictions  are  compared  with  determinations  of  the  ray-path  bending  using 
radiosonde  measurements  of  the  refractivity  profile. 
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SECTION  II 


THE  NBS  SURFACE-CORRECTED  MODEL  ATMOSPHERE 


The  surface-corrected  prediction  techniques  is  based  on  two  empirical 
[3] 

findings  :  on  the  average,  the  change  in  refractivity  over  the  first  kilom¬ 
eter  in  height,  An,  is  related  to  the  surface  refractivity,  N  ,  as  follows, 


-AN  =  7.32  exp  (0.005577  Ng)  ; 


(1) 


and,  on  the  average,  the  radio  refractivity  N  is  an  exponential  function  of 

[3] 

height.  The  1959  NBS  exponential  reference  atmosphere  '  is  defined  by 
Equation  (2): 


N  =  N  exp 
s 


(h-h  ) In 

s 


N 

s 

n  +  an 

s 


(2) 


where 


and  AN 


h  =  altitude  above  mean  sea  level; 

h  =  surface  elevation: 
s 

is  given  by  Equation  (1). 


Figure  1  shows  the  exponential  model  profile  generated  by  a  surface 
refractivity  value  of  354  N  units.  Two  radiosonde  soundings,  both  having 
surface  refractivity  of  354  N  units,  are  also  shown.  These  were  taken  on 
July  21,  1963,  at  Albany,  one  at  1100  and  the  other  at  2300  Greenwich  mean 
time.  It  is  apparent  that  the  actual  refractivity  profiles  differ  from  the  pre¬ 
dicted  profile,  although  the  surface  refractivity  values  are  the  same.  These 
differences  in  profile  mean  that  the  ray-path  bending,  which  is  calculated 
from  the  predicted  (exponential)  profile,  will  differ  from  the  actual  bending. 
The  following  analysis  was  made  to  estimate  the  practical  effect  of  these 
differences  on  height  errors. 


> 
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Figure  1  .  Comparison  of  Radiosonde  N  Profile  with  NBS  Model  Profile 
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SECTION  m 


METHOD  OF  ANALYSIS 


Albany  and  Nantucket  were  selected  for  investigation  because  they  repre¬ 
sent  two  different  climates.  Albany  is  well  inland,  while  Nantucket  experi¬ 
ences  maritime  conditions.  In  fact,  anomalous  propagation  occurs 
extensively  at  Nantucket  during  the  summer  months.  Two  radiosonde 
soundings  were  made  each  day  at  both  sites  with  launch  times  at  1100  and 
2300  Greenwich  mean  time.  Using  the  data  from  these  soundings  in  a  ray¬ 
tracing  computer  program,  the  ray-path  bending  was  then  computed  for  each 
site.  In  the  first  part  of  this  analysis,  the  initial  elevation  angle  of  the  ray 
was  taken  to  be  1  milliradian,  which  represents  a  near-horizon  propagation 
condition.  Under  normal  propagation  conditions,  the  ray  would  be  about 
6500  feet  above  the  ekrth' s  surface  after  traveling  a  distance  of  92  nautical 
miles.  The  Albany  soundings  were  then  used  to  calculate  the  bending  for  a 
ray  with  an  initial  elevation  angle  of  10  milliradians  (corresponding  to  a  ray 
height  of  16,500  feet  at  a  distance  of  124  nautical  miles  under  normal  propa¬ 
gation  conditions). 

For  each  radiosonde  sounding,  the  surface  index  of  refraction  was  used 
with  the  1959  NBS  model  to  calculate  the  bending.  By  comparing  the  data, 
it  was  possible  to  evaluate  the  accuracy  of  this  model. 


5 


SECTION  IV 


RESULTS  OF  DATA  ANALYSIS 

CORRELATION  OF  RADIOSONDE  AND  MODEL  DATA 

The  first  analysis  of  data  considers  the  low-angle  propagation  cases 
where  the  initial  elevation  angle  0 ^  is  1  milliradian.  Figure  2  shows  the 
predicted  and  the  calculated  bending  r  plotted  against  days  during  July  at 
Albany.  There  are  obvious  differences  between  the  curves;  however,  the 
average  of  the  predicted  values  appears  to  follow  the  trend  towards 
slightly  increased  bending  as  the  days  progress. 

A  computer  program  was  used  to  determine  the  auto-correlation 
function  for  the  bending  values  calculated  from  the  radiosonde  profiles. 
This  function  is  shown  in  Figure  3,  and  it  is  apparent  that  the  correlation 
length  is  about  one-half  day.  This  observation  is  somewhat  dependent  upon 
the  one-half-day  sampling  rate. 

The  normalized  cross-correlation  function  was  determined  for  the 
radiosonde  and  model  values  shown  on  Figure  2.  This  function  (seen  in 
Figure  4)  indicates  there  is  a  70  percent  correlation  at  zero  time  delay. 
This  high  degree  of  correlation  indicates  that  the  model  values  closely 
follow  the  variations  in  propagation  conditions. 

EQUIVALENT  HEIGHT-ERROR  DISTRIBUTION 

From  the  large  number  of  ray-tracing  calculations,  a  graphical  re¬ 
lationship  was  established  between  the  overall  bending  r  and  the  elevation 
error  £  .  In  the  range  of  interest  from  about  4  to  14  milliradians  of 
accumulated  bending  (see  Figure  2)  ,  the  relationship  between  r  and  €  is 
approximated  by  a  straight  line  as  shown  in  Figure  5. 
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Figure  2.  Comparison  of  Radiosonde  and  Model  Bending,  Albany,  July  1963,  9  =  1  mr,  Range  of  92  nm 
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Figure  3.  Auto-correlation  of  Radiosonde  Bending,  Albany,  July  1963,  9  =  1  mr.  Range  of  92  nm 
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Figure  4.  Cross-correlation  between  Radiosonde  and  Model  Bending,  Albany,  July  1963,  6  =  1  mr. 
Range  of  92  nm  ° 
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Thus,  the  uncorrected  bending  At  can  be  related  to  an  uncorrected  ele¬ 
vation  angle  error  Ae.  The  height  error  can  then  be  determined  by  multiply¬ 
ing  Ae  by  the  range.  Using  the  above  relationship,  the  uncorrected  bending 
errors  At  between  radiosonde  and  model  values  were  translated  into  a 
height-error  distribution.  Figure  6  shows  this  error  distribution  curve  for 
Albany  during  July  1963.  The  percentage  of  the  height  errors  which  were 
less  than  a  specified  error  I  Ah  |  are  plotted  against  |  Ah  |  .  For  example, 

50  percent  of  all  the  residual  height  errors  are  less  than  +400  feet;  or,  in 
operational  terms,  there  is  50  percent  confidence  that  in  this  situation  the 
height  error  is  less  than  +400  feet.  If  the  confidence  requirement  is  set 
higher,  the  probable  height  error  is  correspondingly  greater. 

ADDITIONAL  ANALYSIS  OF  LOW-ANGLE  DATA 

Similar  analyses  were  made  of  the  low-angle  data  for  Albany  during 
November  and  for  Nantucket  during  July  and  November  1963.  Figure  7 
shows  the  radiosonde  and  model  bending  for  Albany  during  November,  and 
Figure  8  shows  the  cross-correlation  function  for  these  data.  It  is  apparent 
that  there  is  little,  if  any,  correlation  during  this  month.  However,  the 
residual  height-error  distribution  (Figure  9)  indicates  that  the  errors 
averaged  over  November  do  not  differ  essentially  from  those  at  Albany 
during  July  (Figure  6). 

The  comparisons  of  radiosonde  and  model  bending  at  Nantucket  for 
July  and  November  are  shown  in  Figures  10  and  11,  respectively.  The 
radiosonde  results  show  that  anomalous  propagation  is  very  prevalent 
during  July.  The  large  values  of  bending  and  the  frequent  occurrences  of 
trapping  are  evidence  of  the  strong  influence  the  ocean  has  on  such  coastal 
sites.  In  this  type  of  situation,  the  model  which  uses  only  the  surface  index 
is  insufficient  for  making  accurate  calculations  of  bending.  The  height-error 
distribution  curves  for  Nantucket  in  July  and  November  are  shown  in  Figures 
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Figure  6.  Height-Error  Distribution,  Albany,  July  1963,  0  =  1  mr,  Range  of  92  nm 
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Figure  7.  Comparison  of  Radiosonde  and  Model  Bending,  Albany,  November  1963,  0  =  1  mr, 
Range  of  92  nm 
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Figure  8.  Cross-correlation  between  Radiosonde  and  Model  Values  of  r  ,  Albany  November  1963 
9  =  1  mr,  Range  of  92  nm 
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Figure  9.  Height-Error  Distribution,  Albany,  November  1963,  9  =  1  mr.  Range  of  92  nm 
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Figure  10.  Comparison  of  Radiosonde  and  Model  Bending,  Nantucket,  July  1963,  6  =  1  mr, 
Range  of  92  nm  ° 
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Figure  11.  Comparison  of  Radiosonde  and  Model  Bending,  Nantucket,  November  1963,  9  =  1  mr. 
Range  of  92  nm  ° 


12  and  13.  As  Figure  12  shows,  there  is  a  marked  reduction  in  height 
accuracy  due  to  the  anomalous  propagation  conditions  during  July. 

RESULTS  USING  A  HIGHER-ANGLE  CASE 

A  comparison  of  radiosonde  and  model  bending  was  made  using  the 
Albany  data  for  July  but  with  an  initial  elevation  angle  of  10  milliradians 
(Figure  14).  It  is  apparent  that  the  model  follows  the  trend  towards 
higher  bending  as  the  days  progress.  The  corresponding  height-error 
distribution  is  shown  in  Figure  15.  A  comparison  with  Figures  6,  9,  and 

13  shows  that  the  errors  for  0  =10  milliradians  are  about  one-half  the 

o 

errors  with  0  =1  milliradian  and  for  normal  propagation  conditions. 
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Figure  12.  Height-Error  Distribution,  Nantucket,  July  1963,  9  =  1  mr,  Range  of  92  nm 
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Figure  15.  Height-Error  Distributio'n,  Albany,  July  1963,  9  =  10  mr.  Range  of  124  nm 


SECTION  V 


DISCUSSION  OF  RESULTS 


It  is  apparent  from  the  preceding  analysis  that  the  surface-corrected 
exponential  model  provides  a  satisfactory  basis  for  making  first-order 
height-error  corrections  except  under  abnormal  refraction  conditions. 
Although  it  does  not  completely  remove  height  errors,  it  can  be  used  to 
offset  errors  due  to  slow  seasonal  changes  in  bending  conditions. 

Recent  discussions  with  Bean  and  Sweezy  (August  13,  1965,  at  Central 
Radio  Propagation  Lab. ,  NBS,  Boulder,  Colorado)  were  concerned  with  im¬ 
proved  methods  to  obtain  height-error  reductions  at  the  low  elevation  angles. 

The  technique  using  surface  index  and  initial  gradient  data  was  first  reported 
[4] 

in  1963.  As  seen  in  Table  I,  the  attempt  to  extrapolate  this  method  to  the 

low  angles  would  appear  to  provide  results  essentially  similar  to  those 

obtained  by  The  MITRE  Corporation.  However,  Bean  and  Sweezy  stressed 

that  such  an  extrapolation  could  become  erroneous,  and  they  recommended 

reference  to  their  later  methods  for  low-angle  bending  corrections  reported 
[5 1 

in  October  1963. 
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SECTION  VI 


A  COMPARISON  OF  THE  RESULTS  OF  THE  MITRE  CORPORATION 
WITH  THE  NBS  ESTIMATES 


Table  I  shows  the  height  errors  obtained  at  Albany  and  Nantucket  using 

[3] 

the  exponential  reference  atmosphere.  The  68  percentage  level  on  the 
error  distribution  curves  was  selected.  The  listed  height  errors  then 
represent  observations  with  +  1  standard  deviation  of  the  mean  height  error. 
The  extrapolated  results  using  the  NBS  curves  are  also  shown  for  compari¬ 
son,  although  it  is  recognized  that  the  NBS  analysis  was  not  intended  to 

[4  ] 

cover  these  low-angle  cases.  Finally,  the  latest  NBS  results  are  shown 

[5  ] 

using  the  low-angle  correction  techniques.  '  The  locations  analyzed  in 
this  latter  report  were  selected  because  they  represent  climatic  extremes. 

From  the  tabulated  results  of  NBS  and  the  analysis  of  the  13  climatically 

different  sites,  the  use  of  initial  gradient  data  appears  to  reduce  the 

height  error  from  that  obtained  using  the  surface  value,  N  ,  alone.  With  the 

[5]  s 

latest  correction  technique,  the  low-angle  NBS  data  shows  a  marked 
reduction  in  height  error.  The  use  of  initial  gradient  data,  G^,  does  not 
always  contribute  to  an  increase  in  accuracy. 
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Table  I 
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SECTION  VII 


OPERATIONAL  CONSIDERATIONS 


It  has  been  demonstrated  that  by  the  use  of  the  NBS  exponential  refer¬ 
ence  atmosphere,  the  magnitude  and  direction  of  the  average  changes  in 
bending  over  a  one-month  period  can  be  predicted.  The  daily  variations 
were  not  always  predicted;  however,  this  result  was  not  unexpected  due 

[6 

to  the  extreme  variations  which  occur  near  the  earth' s  surface.  Ringwalt 
reported  that  surface  index  changes  as  much  as  15  N  units  in  20  miles  dis¬ 
tance  were  measured  in  Florida. 

Discussions  with  Sweezy  (August  13,  1964,  at  CRPL,  NBS)  indicated 
that  tower  measurements  might  be  the  most  effective  way  to  provide  an 
initial  gradient  measurement  and  a  meaningful  estimate  of  the  surface  index. 
For  example,  a  two-point  measurement  above  the  surface  could  be  extrapo¬ 
lated  to  give  a  representative  surface  value.  These  measurements  might  be 
made  using  the  radar  structure  to  obtain  some  initial  height  advantage. 
Although  this  question  of  minimizing  local  ground  variations  does  not  appear 
to  represent  any  technical  difficulty,  there  is  possibly  an  economical 
problem.  Having  obtained  the  correct  measurement,  there  is  the  remaining 
and  most  difficult  problem  of  how  to  use  the  information  with  available 
field  equipment.  To  apply  the  corrections  would  require  modifications  to 
display  systems  and/or  data-processing  facilities,  which  might  be  central¬ 
ized  so  as  to  serve  several  sites.  These  economical  questions  cannot  be 
answered  without  further  study  and  are  beyond  the  scope  of  this  report. 
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SECTION  VIII 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  results  of  this  investigation,  coupled  with  the  development  of  more 
exact  analysis,  indicate  that  height-error  corrections  methods  are  available 
to  meet  the  requirements  of  most  users.  Excluding  the  need  for  fine 
structure,  for  example,  at  the  Eastern  Test  Range,  the  use  of  realistic 
surface  and  initial  gradient  data  could  reduce  height  errors  to  a  satisfactory 
level. 

There  are  remaining  questions  which  could  best  be  answered  by  an 
error  study  at  a  few  selected  sites.  Such  a  study  would  require  representa¬ 
tive  surface  and  initial  gradient  data  supplemented  by  local  radiosonde  data 
and  aircraft  measurements.  The  extent  of  horizontal  surface  variations 
would  be  investigated  to  determine  if  the  index  measurements  could  be 
meaningfully  applied  in  the  correction  system.  An  aircraft  at  a  known 
height  could  be  used  to  obtain  radar  height  measurements.  The  correction 
technique  would  then  be  applied  to  each  test  mission  to  evaluate  the  method 
and  usefulness  of  the  surface  and  initial  gradient  measurements  in  reduc¬ 
ing  height  errors. 
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